Quantum phases of dipolar spinor condensates 
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We study the zero-temperature ground state structure of a spin-1 condensate with magnetic 
dipole-dipole interactions. We show that the dipolar interactions break the rotational symmetry of 
the Hamiltonian and induce new quantum phases. Different phases can be reached by tuning the 
effective strength of the dipolar interactions via modifying the trapping geometry. The experimental 
feasibility of detecting these phases is investigated. The spin-mixing dynamics is also studied. 
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A dipolar condensate is a condensate possessing the 
dipole-dipole interactions, in addition to the usual s- 
wave contact interaction. Recent theoretical studies have 
found that, the long-range dipolar interactions (which 
may originate either from the intrinsic atomic Q] or 
molecular 0, 01 Q| magnetic dipole moments or from the 
induced electric dipole moments la) may greatly affect 
the stability of the condensate laMEll jl an d significantly 
modify the excitation spectrum |& lid. llH Il2j. Further- 
more, the dipolar interactions are responsible for the gen- 
eration of translational entanglement EH an d a variety of 
novel quantum phases [Tlj in condensates confined in op- 
tical lattices. All these works have focused on scalar con- 
densates. The effect of the dipolar interactions on spinor 
condensates confined in multi-well potentials is explored 
in the works of Refs. Eil> where only the the dipolar in- 
teractions between different potential wells are taken into 
account. In this Letter, we study the properties of dipo- 
lar spinor condensates in a single trap. We show that the 
interplay between the spin-exchange and the dipolar in- 
teractions gives rise to extremely rich physics unexplored 
so far and that as the effective dipolar strength depends 
on the condensate aspect ratio, transitions between new 
quantum phases resulting from the dipolar interactions 
can be induced by modifying the trap geometry. 

The intrinsic magnetic dipole moment of an atom is re- 
lated to its total angular momentum as fip = — pj?//flF, 
with gp being the Lande g-factor, /is the Bohr magne- 
ton, and F the angular momentum operator. Thus the 
dipole-dipole interactions will affect the condensate once 
the spin degrees of freedom become accessible as in an 
optical trap. Since the first creation of spinor conden- 
sates in dilute alkali atomic vapors |l6[ . we have wit- 
nessed tremendous experimental and theoretical studies 
on these systems [H El E 0, HH IH H . The key 
feature of the spinor condensate is the spin-exchange in- 
teraction which is responsible for a variety of interesting 
physics such as the spin-mixing |21|. phase conjugation 
[23|. spin domain formation and topological defects 
generation |24|. All these studies, however, neglected the 
dipolar interactions, except in Ref. |25j, where Gu sug- 
gested that the internal dipolar field may induce sponta- 
neous magnetization in a spinor condensate. Neverthe- 



less, the dipolar interactions are not treated explicitly in 
that work. 

Consider a spin-I condensate with TV atoms. The 
Hamiltonian without the dipolar intereactions, in second 
quantized form, reads [2(j 
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where M is the mass of the atom, i> a ( r ) (a = 0, ±1) 
denotes the annihilation operator for the mp — a com- 
ponent of a spin-1 field. The trapping potential V ey ± is 
assumed to be spin-independent. The collisional inter- 
action parameters are Co = Anh 2 {a,Q + 2a2)/(3M) and 
c 2 = 4:nh 2 (a 2 - o )/(3M) with a f (f = 0, 2) being the s- 
wave scattering length for spin-1 atoms in the combined 
symmetric channel of total spin /. For the two exper- 
imentally realized spinor condensate systems ( 23 Na and 
87 Rb), we have |ca| Ceo. 

The Hamiltonian for the dipolar interactions reads 
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where Cd = ii<iy? B g 2 p I (Ait) and e = (r— r')/|r— r'| is a unit 
vector. The total Hamiltonian is then -f/tot = ^sp + -ffdd, 
from which we can derive the Hesenberg equations for 
the boson field operators. In the standard mean-field 
treatment, one replaces the field operators ip a by their 
expectation values <f> a = (i[) a ), which yields the so-called 
Gross-Pitaevskii equations (GPEs). 

Before proceeding further, let us first estimate the rel- 
ative strengths of the spin-exchange and the dipolar in- 
teractions which are characterized by C2 and Cd, respec- 
tively. If we take eto = 50as (as being the Bohr radius) 
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and a 2 = 55a B for 23 Na [H, we find c d /\c 2 \ ^ 0.007. 
Similarly, for 87 Rb (a = 101. 8a B and a 2 = 100.4a B 
|22^). we have Cd/|(^| — 0.1. Hence at least for 87 Rb, 
the dipolar interaction is not negligible compared to the 
spin-exchange interaction. 

The facts that \c 2 \ <C Co and Cd -C | C2 1 ins pire us to in- 
voke the single mode approximation (SMA) |2]j, namely, 
V'aM ~ </ , (r)a Q , where <^>(r) is the spin-independent con- 
densate spatial wave function, a a is the annihilation op- 
erator for mp = a component. The validity of the SMA 
can be checked by solving the GPEs numerically for the 
ground state wave functions <j> a . We have performed 
these calculations and found that the SMA is valid for 
Cd < . 2 1 C2 with scattering lengths of 23 Na and 87 Rb. 

The Hamiltonian H sp under the SMA is given by (after 
dropping spin- independent constant terms) plj : 
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(1) 



where c' 2 = (C2/2) J cZr|</)(r)| 4 is the spin-exchange inter- 
action coefficient, and L = atF^a^ is the total many- 
body angular momentum operator of the system, whose 
eigenstates and eigenvalues are defined by 

L 2 |Z, m) = 1(1 + 1)|Z, to), L z \l, to) = m\l, to), 

where m = 0, ±1, . . . , ±Z and for a given total number of 
atoms N, the allowable values of I are I — 0, 2, 4, . . . , TV 
for even N and I — 1, 3, 5, . . . , N for odd N. The ground 
state properties of the spinor condensate under JJ sn i s 
therefore completely determined by the sign of c' 2 [21j : 
For c' 2 > (antiferromagnetic), the ground state is given 
by the spin singlet \G) = |0,0); for c' 2 < (ferromag- 
netic), it is given by \G) = \N, m) which has a (2 AT + 1)- 
fold degeneracy with to taking any integer numbers be- 
tween —N and N. 

The dipolar interaction Hamiltonian under the SMA is 
given by 
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In general, this is still in a very complicated form. Re- 
markably, for a condensate with axial symmetry (which 
happens to be the most experimentally relevant case), 
with its symmetry axis chosen to be along the quantiza- 
tion axis, z, H^d takes a very simple form: 



H, 



dd 



-^L 2 + 3c^L 2 + n ), 



(3) 



where no = o,q&o is the number operator for mp = 
and c' d = (cj/4) / / drdr' 1 0(r ) | 2 1 0(r' ) | 2 ( 1 - 3 cos 2 e )/ |r- 
r'| 3 with 9 e being the polar angle of (r — r'). The total 



Hamiltonian under the SMA is then 

fftot = (4 - c' d )t 2 + 3c' d (L 2 z 



n ) 



(4) 



Before we discuss the ground state of -fftot , we want to 
point out that H sp in possesses SO(3) symmetry, i.e., 



it is rotationally invariant in spin space |20j. The pres- 
ence of the dipolar interaction, however, breaks this sym- 
metry. This is not unlike the situation in superfluid 3 Hc, 
where the dipolar interaction between nuclear spins, de- 
spite of its smallness, breaks the spin-orbit symmetry and 
is crucial for the understanding of the superfluid phases 
of the system 26]. As we shall see below, the dipolar 
interaction of atomic spins also results in new quantum 
phases in spinor condensate. 

From Hamiltonian |J2J| , one can expect that the behav- 
ior of the dipolar spinor condensate should be very sensi- 
tive to the signs of both (c' 2 — c' d ) and c' d . Unlike c' 2 , both 
the sign and the magnitude of the dipolar interaction 
coefficient c' d depend on the geometric shape of the con- 
densate, which makes the system highly tunable through 
the modification of the trap aspect ratio. As an example, 
we assume the condensate wave function has a Gaussian 



form 4>(r) — n 3 / 4 k 1 / 2 
calculation shows that 



c'M\ = 27rc dX (K)/(3\c 2 \), 



, simple 
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where x( K ) is a monotonically increasing function of the 
condensate aspect ratio k, bounded between —1 and 2, 
and passing through zero at k = 1 [27j • From this it is 
clear that both the sign and the strength of the effective 
dipolar interaction can be tuned with trapping geometry. 

To gain more insights into the structure of the ground 
state, let us first neglect the no-term in -f/tot ■ The remain- 
ing Hamiltonian, denoted by H [= (c' 2 — c' d )L 2 + 3c' d L 2 ], 
has a diagonalized form under the basis states \l,m). The 
ground state of Hq can be easily found and the phase di- 
agram is illustrated in Fig. 2] We see that the ground 
state of Hq is divided into three regimes in the c' 2 -c' d pa- 
rameter space, and the ground state angular momentum 
quantum number I can be either N or 0. 
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FIG. 1: The ground states of H , where l = N if N > 
(c— l)/(2c+l), otherwise lo = 0. For typical parameters, the 
inequality is usually satisfied, so we can always take lo = N. 

For the case of I = N, we expect the ground state of 
the total Hamiltonian H to t to be very close to that of 
Ho, since the no-term is at least a factor oil/N smaller 
than the rest. In contrast, for the / = case, the no-term 
is expected to be important and the true ground state 
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should deviate from |0,0) significantly unless c = 0. We 
found the ground state of H tot numerically by expanding 
the Hamiltonian onto a Fock state basis \N-%, No, N-\) 
with the constraint N — ^a, an( ^ the calculations 
confirmed our qualitative arguments above. In the fol- 
lowing, we will discuss the properties of the ground state 
in more detail and we will present our results for the cases 
c' 2 > and c' 2 < separately. 

c' 2 > case: Figure |3 (a) illustrates the overlap be- 
tween the true ground state (denoted by \G)) and that 
of Hq ■ The overlap is unit in the regions c < — 1/2 and it 
quickly approaches unit for c > 1. For —1/2 < c < 1, the 
overlap decreases rapidly from 1 when c deviates from 
0. In this region, in general \G) = ^;<7i|J,0), with the 
coefficients gi dependent upon c, is a superposition of 
different angular momentum states with (L z ) = 0. The 
normalized ground state population and its variance are 
presented in Fig. [21(b) and (c). Two sharp phase bound- 
aries can be seen at c = —1/2 and c — 1 in agreement 
with Fig.^ The variance is presented in terms of Mandel 
Q-factor. Q < 0, = 0, and > represent sub-Poissonian, 
Poissonian, and super-Poissonian distributions, respec- 
tively. In particular, Q = — 1 represents a Fock state 
with vanishing number fluctuations. 

For c < -1/2, \G) = \N,±N) is a Fock state with 
all the population in either rap — 1 or —1 state. This 
is the state with spontaneous magnetization discussed in 
Ref. [25| . Since the ground state in this regime is two- 
fold degenerated with opposite magnetization, a chain 
of such condensates will form an example of the Ising 
model of statistical mechanics. For c > 1, \G) ~ |JV, 0). 
Expanded onto the Fock state basis \N—i, Nq, N\), we 

have \N,0) = ^k=o 9k\N/2-k, 2k, N/2-k) (assuming TV 
to be even), where for N ^> 1, the expansion coefficients 
g k = (8/A^7r 2 ) 1 / 4 exp[-4(/c - N/A) 2 /N], from which we 
immediately have (no) = 2(n±i) = N/2 , Qo = —1/2 and 
Q±i — —3/4. Hence the populations in all three spin 
components have sub-Poissonian distributions. Finally, 
for —1/2 < c < 1, population in mj? = state has a 
super-Poissonian fluctuation with maximum fluctuations 
occur at c = 0, while the population fluctuations inmj? = 
±1 approach the Fock state limit for positive c. 

c' 2 < case: For this case, the overlap between the true 
ground state and that of H is nearly perfect. The phase 
diagram is characterized by a sharp boundary at c = 0, 
and the ground state is given by the maximally polar- 
ized state \N, ±N) or \N, 0) for c < and c > 0, respec- 
tively. Since the phase boundary occurs at c = (vanish- 
ing dipolar interactions), this phase transition should be 
readily verifiable in Rb spinor condensates by changing 
the condensate shape from cigar to pancake, or vice versa. 
As we have mentioned earlier, without the dipolar inter- 
action, the ground state of a ferromagnetic spin-1 conden- 
sate is given by \N, m) for any integer m between — N and 
N (i.e., the spin vector has no preferred spatial orienta- 
tion) . Now the dipolar interaction breaks this degeneracy 
and orients the spin vector along (perpendicular to) the 
axial direction for a cigar-shaped (pancake-shaped) con- 
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FIG. 2: The c-dependence of the overlap between the true 
ground state \G) and that of Ho (a), the normalized ground 
state population (b), and the number fluctuation (c). The 
solid and dashed lines in (b) and (c) denote the rriF — 
and —1 components respectively. The inset in (c) shows the 
detailed behavior of Q a at the vicinity of Q a = 0. The results 
presented here are for c' 2 > and N = 100 l2Sll . 



densate, a clear manifestation of the symmetry-breaking 
properties of the dipole force. 

In addition to the ground state phase structures, the 
spin-mixing dynamics |2l| of the dipolar spinor conden- 
sate can also be studied by numerically evolving an initial 
state under H tot . Figure [3] illustrates one example. Here 
the initial state is given by the Fock state \N/2,0,N/2) 
with half population in mp = 1 and —1 components, re- 
spectively. The spin-mixing dynamics will quickly drive 
the system into a quasi-steady state. Again we want 
to focus on the effect of the dipolar interactions whose 
strength depends on the condensate aspect ratio k. In 
Fig. the steady state population of spin-0 is plotted as 
a function of c, while the inset replots the population as 
a function of k using parameters for 87 Rb. As we can 
see, modifying the condensate aspect ratio changes sig- 
nificantly the population distribution in the steady state. 

In summary, we have studied the ground state proper- 
ties and spin-mixing dynamics of a dipolar spinor conden- 
sate. We have shown that the dipolar interaction breaks 
the rotational symmetry of the Hamiltonian and as a re- 
sult the ground state is characterized by several distinct 
quantum phases depending on the relative strengths of 
the spin-exchange and dipolar interactions. The transi- 
tions between these phases may be induced by simply 
modifying the trapping geometry of the condensate. 

We have neglected any external magnetic field in our 
study here. The presence of external fields will affect the 
orientation of the spin and hence change the phase di- 
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FIG. 3: The c-dependence of the steady state population of 
rriF = component for the initial state \N/2,0, N/2) with 
TV = 10 4 and c'2 < 0. The inset indicates its jt-dependence for 
87 Rb condensate with a Gaussian wave function. 



agram. For the change to be insignificant, we need to 
reduce the field strength such that the zeeman energy 
is weaker than the dipolar energy. For typical values of 
alkali atoms, this requires to control the magnetic field 
below 1CT 4 Gauss |2<|. This requirement will pose an ex- 
perimental challenge (none of the experiments on spinor 
BEC so far 0,|l|BlI3 has met this field-free re- 
quirement), but is within reach with current technolo- 
gies. We hope our work will stimulate experimental ef- 
forts along this line. The effect of the magnetic field on 
dipolar spinor condensates is currently under study. Our 
future works will also include the study of spin-2 dipolar 
spinor condensate 0, 0] . We believe that these stud- 
ies will open many unexplored and promising avenues of 
research in the field of quantum degenerate atomic gases. 

We thank Jason Ho and Michael Chapman for helpful 
discussions. 
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